Entropy of gravitating systems: scaling laws versus radial profiles 
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Through the consideration of spherically symmetric gravitating systems consisting of perfect fluids 
with linear equation of state constrained to be in a finite volume, an account is given of the properties 
of entropy at conditions in which it is no longer an extensive quantity (it does not scale with system's 
size). To accomplish this, the methods introduced by Oppenheim [1] to characterize non-extensivity 
are used, suitably generalized to the case of gravitating systems subject to an external pressure. In 
particular when, far from the system's Schwarzschild limit, both area scaling for conventional entropy 
and inverse radius law for the temperature set in (i.e. the same properties of the corresponding black 
hole thermodynamical quantities), the entropy profile is found to behave like 1/r, being r the area 
radius inside the system. In such circumstances thus entropy heavily resides in internal layers, in 
opposition to what happens when area scaling is gained while approaching the Schwarzschild mass, 
in which case conventional entropy lies at the surface of the system. The information content of 
these systems, even if it globally scales like the area, is then stored in the whole volume, instead of 
packed on the boundary. 
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I. INTRODUCTION 



The study of the connection between entropy and gravity, performed even simply at a classical level, has 
deserved in recent years very interesting results. Examples have been found in which the conventional (as 
opposed to black hole) entropy of gravitating systems scales like their area and the temperature like l/R (being 
R their area radius), that is the same as the celebrated laws of black hole entropy and temperature if area and 
radius refer to horizon 0, Q ■ 

Black hole properties are apparent at a semi-classical level in which quantum fields live in a continuous 
^ 'j curved spacetime, describing gravity; these properties moreover, in particular the area scaling law of black hole 
?H . entropy, are considered to give hints or constraints on a fundamental description of spacetime itself, possibly at 
IPP" a quantum level, of a holografic nature 0, @ . Now, the fact that purely classical systems can exhibit analogous 
^ , scaling laws in the corresponding conventional thermodynamical quantities is intriguing. The study of these 
systems could shed light on what in the area scaling property of black hole entropy truly points to a fundamental 
rS | origin. 

In Q gravitating systems consisting of perfect fluids with linear equation of state p = np, with p and p the 
pressure and energy density and k a (non- negative) constant, have been considered. The authors have found 
the remarkable property that for the stiffest case compatible with causality, namely the k = 1 case (see Q), the 
equilibrium configurations can exhibit area scaling of their entropy, as well l/R scaling of their temperature, 
suggesting a possible connection between holografy and causality. This recalls the formulation of the holographic 
principle, as given by [8[, which has causality built-in. In Q, the entropy of a spherically symmetric gravitating 
system in equilibrium has been shown to become area scaling as the radius approaches the Schwarzschild value; 
in this limit the total entropy of the body is accounted for only by the entropy of the more external spherical 
shells. 

In general when gravity enters the game, extensivity of entropy is lost even if, by the equivalence principle, 
it is mantained locally. In [l| a method has been introduced to study the degree of non-extensivity of entropy 
(its scaling behaviour with system's size) for systems with long range interactions (other approaches to non- 
extensivity, trying to characterize and understand it outside standard thermodynamics, use Tsallis entropy [l(J; 
see also [13). Applying this method to a particular system consisting of a spherically symmetric distribution 
of perfect fluid with constant density, it appears that at increasing M/R (M is the mass, or mass-energy, as 
spelled out in [l2j], known also as ADM mass) entropy starts going towards area scaling, with an increasing 
contribution from external layers (the limit R — 2M is however not reachable as equilibrium is possible only if 
R > (9/4)M (see, for instance , fl^ )). in agreement with what expected from Q. The idea in the present paper 
is to go back to the systems [6|, known to reach conditions with area scaling of entropy, and to investigate, 
through a generalization of the method introduced in [3 , the dependence of non-extensivity on M/R as well to 
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study the evolution of the entropy profile. 

The paper is organized as follows. In section II we examine the equilibrium configurations of the fluids 
under consideration first under conditions in which they are unconstrained and next when they are instead 
constrained to be in a fixed volume. For these latter conditions in section III we try to characterize the degree 
of non-extensivity of the systems following the approach [l[ , trying to extend it to cover the case in which an 
external pressure is applied, i.e. our case. A global temperature is introduced and an expression is found that 
for each k gives the degree of non-extensivity of the entropy of the system for assigned (central) energy density 
and size. In section IV first we use this result to investigate the scaling laws for entropy; next the scaling laws 
for global temperature are also derived and, finally, the entropy radial profiles are considered, in their relation 
with the scaling laws of entropy, in particular for the most interesting case k = 1 . We conclude with a summary 
of our results in section V. 



II. EQUILIBRIUM CONFIGURATIONS FOR LINEAR PERFECT FLUIDS 

As in I &] , let us consider a perfect fluid at thermodynamic equilibrium with linear equation of state 



p = np : (1) 

where p and p are respectively the energy density and pressure in the rest frame of the fluid element; k is a 
constant with < k, < I where, as said, K — 1 corresponds to the stiffest conditions compatible with causality. 
The stress-energy tensor is 



T a /3 = (p + p)u a Uf3 + pg a/3: (2) 

where u a is the fluid element 4-velocity and g a /3 is the metric tensor. In the case of static equilibrium con- 
figurations (necessarily with spherical symmetry), from (fTJ) and ([2]) Einstein field equations in Schwarzschild 
coordinates can be written 
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where the derivative symbol means derivative with respect to radial coordinate (area radius) r and where the 
metric is 



ds 2 = -f(r) dt 2 + h(r) dr 2 + r 2 dfi 2 , (4) 

with / and h two non-negative functions of r, to be determined. 

The properties of the solutions to these equations are known since a long time (l3| . The equilibrium con- 
figurations are completely determined b y th e equation of state (fTJ); they can be derived solving directly the 
Tolman-Oppenheimer-Volkoff equation (jl4l. [l5| and [IH). Aiming however at obtaining entropy profiles, we 
want to stress the dependence on temperature. We choose then to solve equations ((3|) inserting explicitly this 
dependence since start. For static gravitational fields, Tolman relation L l(3] prescribes that T^J — goo be constant 
throughout the system, where T is the temperature in the rest frame of the fluid element. This means for us 
that 

T{r)^W)= const. (5) 
For our fluid, thermodynamics gives the following scaling laws [13, EH 



p cx T (K+1)/K , 
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with s the entropy density. The proportionality constants depend on the particular realization of the fluid. 
Applying these scaling laws locally throughout our system (equivalence principle), from equation (O we have 
thus 

p(r)f^ (r) = const = — , (8) 

with C some constant with the dimensions of energy density (length -2 in geometrized units). Note that while 
the value of C depends on the chosen unit for length, the quantities p/C and r\/~C are dimensionless. 
Making use of ([8]), the first and the second of equations ([3]) become 

C b! h - 1 

(9) 



/t^t rh 2 hr 2 
C f h-1 



(10) 



while the third, as can be verified (and as noticed in [jjj], in the case of a photon gas), is always identically 
satisfied if equations (Pl fTUj) are satisfied, as a confirmation of Tolman result. 

Before we start to solve these coupled differential equations, let us study some general characteristics which 
the solutions f(r) and h(r) should satisfy. h(r) can be put in the form 



h(r) = 1 ; w , (11) 

1 1 l-2m(r)/r' y } 



(see for example |12| ) with m(r) given by 



C r f 2 

m(r) — I 47TT 2 p(r)dr — — / — r^dr, (12) 
' 2 J fB( r ) V ; 



o 

that is, m(r) is the ADM mass inside r. Here we have introduced the function 

= (13) 

so that g(n) is monotonically decreasing in the relevant range < n < 1, with g{n) > 1 and g{n) — > oo for n 
approaching 0. The cases of perfect fluids with k — (dust), k = 1/3 (blackbody radiation) and n = 1 (stiff 
matter) correspond respectively to g = oo, g = 2 and g = 1. In the following we will often refer to g instead of 

AC. 

For each allowed equilibrium configuration, in the limit r —> 

m(r) ~ ^7rr 3 p(0), (14) 
so that m(r)/r — > 0. This implies /i(0) = 1 for each solution with p(0) finite. As regards /, from |8]) we have 

P(0)/ 9 (0) = ^, (15) 

so that finite p(0) implies a non zero /(0). Taking the derivative of (fTTj) and using (fT4|) and (fl~5|) . in the limit 
r -> we have h'(r) ~ x7r/9(0)r and /'(r) ~ r/(0)^7r^_p(0), so that /i'(0) = and /'(0) = 0. 
We want now to solve equations ([9l fT0|) . We can write equation (fTO]) as 

fg-i(f + r f) 

*-<*-» d. + { fc-w (16) 



and, substituting this expression of h in equation @, we obtain 



4 



,„ = 2(g + l)Cr/ 2 + 3gQ 2 //' - 2(2g - l)P +1 f + gCr 3 /' 2 

1 CV3/ + (2g-l)r/9+i 1 ' 

so that equations (|16M17p are equivalent to the starting equations ([51 fTU|) . Using, in addition, the previous results 
about the behaviour of / near the origin, the problem is thus reduced to find the solutions / of the second order 
equation (JTTJ) with initial conditions 



/(o)-(^oj) 1/9 ; /'(o) = o- (is) 

For assigned g, for each choice of the pair C, p(0) we have a solution f(r). Note however from equation 
([8]) that physical configurations p — p(r) depend only on the ratio Cj f 9 and not on C and / separately. This 
implies that for any chosen p(0), the (infinite number of) solutions / that we find when C is varied, correspond 
always to the same physical configuration p — p{r). For each g, to fix p(0) corresponds to uniquely fix the 
physical configuration, irrespective of the separate values of C and /. Without loss of generality we can then 
study the solutions to equation (fT7|) for a single arbitrary choice of the value of C 0); let us put for simplicity 
C = 1. All possible static configurations of our system are then in 1-1 correspondence with all the solutions / 
of equation 



,„ = 2(g + l)r/ 2 + 3gr 2 //' - 2(2g - l)P +1 f + gr 3 f 2 

J r 3 / + (2ff _l) r/9 +l ^ 



with initial conditions 



f^)=( 1 ^ m ) 1 '\ /'(0) = 0, (20) 

where p(O) is any positive number and with h given by (|16[) . Due to ©, each specification of p(0) is equivalent 
to a specification of the central temperature Tq ■ 

Some symmetry considerations can now help to understand how the set of solutions / (and h) is arranged. 
Consider the transformation 



r -> r = \r (21) 
f^f = \ 2 /3f. (22) 

We have 

/(f) = A 2 /*/(f/A) 

or 

/» - X 2/9 f(r/X) 

and then 

f'(r) = Af-W/A) 

and 

/"(r)=Al- 2 /"(r/A) 

so that, if f(r) is a given solution, f(r) is also a solution, as can be easily verified going through equation (fT9|) . 
f(r) is such that /(0) = A 2 ^ 9 /(0). On the other hand we know that, starting from a given solution /i, we 
can obtain all solutions for the assigned g if we change the initial conditions, /i(0) — > f\(0) = A/i(0), with A 
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FIG. 1: / versus r for k = 1, i, § (highest: k = 1, lowest: /t = ~) for initial conditions (|26|l . see text. 



spanning all positive real values. Thus, in terms of an assigned solution fi(r), a function f\ is a solution if and 
only if 



with A > 0. For the solutions h we obtain 



A(r)=A/i(r/A*) 



/ lA (r) = / H (r/A^), 



(23) 



(24) 



where /iaM and ft-i(r) corresponds to f\(r) and through (fl~6)l . We see that, from {§]), the transformations 
we are considering are equivalent to say that equations ((3|) are invariant under the transformation [(| HO, E[ 



r — > Ar 
P - P/A 2 
p^p/A 2 . 



(25) 



Summing up, once we know the behaviour of solutions / and h for a given particular initial value /(0) = a\ 
of /, arbitrarily fixed (> 0), we can readily get the corresponding behaviours of the solutions for every other 
assigned initial value a,2- Let us study for simplicity the case /(0) = 1, that is the solutions to ([i"9]) with initial 
conditions 



/(0) - 1; /'(0) = 0. (26) 

From numerical integration of equation (|19p and from equation (|16p we obtain the solutions / and h reported 
in Figures [T] and [2] for the cases g = 1, 2, 5 (i.e. k = 1, ^, ). The structure of the corresponding equilibrium 
configuration can be seen in Figures [3] and [U that respectively report p and 2m/r, as a function of r, for the 
same chosen values of g. Our choice C = 1 implies that p and r are given here respectively in units of C 
and C -1 / 2 . From (JSJ and (|23|) the properties of p for every other solution can readily be inferred. Note that 
p(r) approaches asymptotically as r — > oo. This implies that for this configuration and, by scaling property 
(|23"1) . for every other solution, the radius of the system is infinite as expected. These radial profiles agree with 
the corresponding results reported in [22| within an extensive study of the conditions of general relativistic 
instability of finite spheres consisting of exactly the same fluids we are examining here. 

Let us consider in some detail the functions h and 2m jr (Figures and 01 . h starts from 1 at r = 0, and 
2m /r from 0; they reach their maxima, after which they go towards their own limiting values [l9l |. 



h(oo) = 2, 



7 34 
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respectively for g = 1, 2, 5, shown in the Figures, with damped oscillations around them, in which both the 
amplitudes of the deviations from the limiting values as well their derivatives at any order reduce asymptotically 
to 0. Analytically both limiting values can easily be verified. In fact consider for example that from (fTTj) we 
have 



^ = iLll. (29) 
r h 

In equation ([9]) we see then that, at extremal points, that is when r has values r* such that h'(r*) = (that is 
(2m(r) / r)' r=r » = 0), we have 

r* ~ f 9 (r*y m 

As extremal values tend to (2m/r) oc when r — > oo, from equation (JTUJ) considered in this limit, using (f!?9")) we 
have 

1 /2m \ 2 ( (2m\ \ (2m 



2q — 1 \ r I q \ \ r J m J \ r 
so that 



2m\ 4-2/g 



and 



4 - 2/ 9 + 2.g 

4 - 2/ 5 + 2.g 



(31) 



ft(r) = .„ . = ^ - (32) 

w 1 _ 2™ 2g V ; 

and thus the values (|28|) , as well (|27|) . Here we made use of the fact that in the limit r — » cxd, (r 2 / f 9 )' — implies 
/' = jjf/T. Note that equations ([2^)1 and ([2^|) determine 2m/r for every solution; we see that 2m/r scales as h. 
From scaling property (|24[) and from equations (|20p and (|22|) . the asymptotic behaviour of h and 2m /r when 
r — > 00 implies that in the limit of very large central density p(0) or very large Tc, h(r) ~ const = h(oo) and 
2m(r)/r ~ const — (^) , for every r, except the region of r very near to 0, in which /i goes from 1 to h(oo) 
and 2m I r goes from to j with damped oscillations. 
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3: p versus r for k = 1, i, | (highest: « = ~, lowest: /t = 1) for initial conditions (|26[) . namely p(0) = 1/8tt, see 



As regards /, from ifTU)) we see that f'(r) > 0, Vr and from ([50)) we have that in the limit r — ► oo, 



Due to equation ([2"3"[) this implies 



and then 



hm / A (r) 
lim A ( 



lim f x {r) 
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T Um^/ A (r) = UmA/(^; 
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(34) 



that is when the central density is allowed to increase without limit, the density at each fixed radius r does not 
increase unlimitedly but approaches a limiting value proportional to r~ 2 . Looking at (|31| . the proportionality 
constant here is such that these limiting conditions correspond to the analytic (unphysical) solution for the 
pressure (or energy density) [H[ 



P 



1 



2tt(1 + 6k + k 2 ) r 2 



(35) 



in agreement with what stated in [a, [l3, 22]. This analytic solution is represented in Figures [5] and [J] by the 
horizontal lines; the intercepts of the physical solutions with the analytic one follow a well known geometrical 
progression with ratio depending on k [13l l22| . 

Finally, from scaling properties (|23H24|) . from initial conditions (|2U)l and from the behaviour of f(r), h{r) 
and p(r) near r = 0, shown in Figures [TJ [5] and [31 we have that in the limit of very low central density, f(r), 
h(r) and p(r) are constant in a ever increasing range of r, so that the flat-space limit with constant densities is 
recovered. 

To obtain now a system with finite size we can imagine to enclose our fluid inside a perfectly insulating 
spherical cavity, with area radius R, and to assume that thermodynamical equilibrium is reached, that is the 
temperature of the walls is equal to the temperature of fluid at r = R. This will correspond to a certain 
central energy density pc for the cavity. Due to the local uniqueness of solutions to equation (|19|) with initial 
conditions ([2"Uj) . this is equivalent to take a certain (physical) solution to equation (fTO)) (that with p(0) = pc) 
and to cancel all fluid outside r = R. As regards energy density radial profile p(r), this means that there is a 1-1 
correspondence between all admitted cavities filled with linear fluids (that is every admitted central densities as 
well radii) and all possible p(r) distributions in finite ranges [0,R], where p(r) = l/(8irf 9 (r)) and / is solution 
to p9[) with initial conditions ([2"0)l . 
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FIG. 4: 2m jr versus r for k = 1, |, | (highest: « = 1, lowest: n = i) for the same conditions as previous Figures. 
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FIG. 5: 2M/R versus z = p(0)R 2 (R can be thought fixed), for k = 1, i, | (highest: k = 1, lowest: k = |" 



As allowed by Birkhoff theorem [TJ, |24j , if the constraining system is assumed spherically symmetric and 
with finite extension, the Schwarzschild coordinates can be chosen in such a way that outside a certain radius 
R* exactly the Schwarzschild form of the metric obtains 



ds 2 = -(1 - 2M*/r) dt 2 + dr 2 + r 2 d£l 2 . (36) 

v 1 ' l-2M*/r y ' 

The parameter M* here is precisely the mass-energy or ADM mass for the whole system, that is fluid + 
constraining system. Due to our symmetry conditions we can look at the quantity 

M = ^r 2 p(r)dr = - j J^ dr > ( 37 ) 

as the mass due to the fluid alone, that is the mass would be probed by Kepler-like experiments with gravitating 
bodies far away from the fluid alone. It represents a meaningful quantity, proper of the ball of fluid contained 
in the cavity, independent of the characteristics (mass, extension,..) of the constraining system. As the metric 
outside has been explicitly chosen, the continuity of the metric tensor components g a p everywhere, even across 
boundaries, fixes, for every assigned pc and g, f(r) Vr. Due to (TTS]) the value of C is fixed also. In other 
words, once a single coordinate system is explicitly chosen, the arbitrariness in the choice of the constant C in 
equations (fi^- (fT0)) is lost. 

Fig. [5] permits to see, for a cavity with fixed area radius R at equilibrium, the dependence of its mass on the 
central density or temperature, for different assigned g. It reports 2M/R as a function of z = p(0)R 2 . The plot 
is done in terms of z instead of p(0) = y in order to obtain for 2M/R an universal shape, i.e. not depending 
on the chosen R. The expression for z indicates how to rescale the abscissa to obtain, from the given plot, 
the actual relation between 2M/R and y for each fixed value of R. The explicit form of the plotted function, 
expressed in terms of the solution hi for initial conditions (l^rj)) (the solution h presented in Fig. [5]), is 
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2M(z) = 

as can be easily derived from equation (12"9")) and from the scaling property (12"4")) . A maximum (2M / R) max = 0.544 
is found (with the extremal point at z extr = 0.655), corresponding to the case g = 1 (k = 1); it gives the 
maximum fraction of the Schwarzschild mass a perfect linear fluid at equilibrium confined inside a given radius 
can have. The maximum value for the case g — 2 (k = 1/3), namely (2M/ R) max = 0.493, is in agreement 
with [20] . All these results agree with what reported in [13] , where in addition a prescription is given to relate 
(2M/ R)max values when k is varied. 

In Fig. [5] we see that once the extremal central energy density or temperature are overcome, equilibrium 
configurations are still possible, at lower values of M. We have thus larger central temperatures, which give rise 
however to a lower total mass for the system. This is by no way a problem for energy conservation, obviously. 
The total energy of the system, given by M, contains in fact both the component of proper energy p and the 
(negative) component of gravitational potential energy [T2I |25| . Evidently, when the extremal temperature is 
overcome, the system admits a new equilibrium configuration with a stronger gravitational energy component. 

A detailed discussion of the stability of the confined fluids we are considering here is in [22]. According to 
those results a new mode on instability occurs whenever a (non trivial) local extremum in our Fig. [5] is present. 
This amounts to say that, as far the equilibrium of the fluid alone is considered (neglecting the constraining 
system [Iff]), starting from a situation of stability at low z (very low energy density and/or very small radius), 
the first maximum marks the transition to instability and this latter will never desappear for larger values of z, 
so that stability is for z < z extr . This agrees with the results in [2C| corresponding to the case k = | and with 
the general discussion of thermodynamic equilibrium and stability conditions in [27| • 



III. NON-EXTENSIVITY 



The aim of this section is to characterize the systems under study in terms of their degree of non-extensivity. 
This is done using the methods introduced in generalizing them to the case of gravitating systems subject 
to an external pressure. 

The first law of thermodynamics as applied to an infinitesimal volume dV of a homogeneus system reads 

dE = T dS -p dV + [i dN, (39) 

where dE, dS and dN are respectively energy, entropy and number of particles in the volume dV, and \x the 
chemical potential. Integrating this equation to a tiny volume V means to sum over elementary volumes, 
each possessing always the same values of the intensive quantities. If the quantities on which we integrate are 
extensive, that is if they scale with size, integration gives 



E = TS~pV + pN, (40) 

so that this relation (known as Gibbs-Duhem relation) requires and manifests extensivity. This can conveniently 
be expressed as 1 = or 1 = or 1 = E+ ^f_^ N if, respectively, no p and p terms, p terms or p and 

fi terms are chosen to be included in the expression of S. In particular, without p and fj, terms in S, a linear 
relation 



S(XE) = XS(E) (41) 

is implied, where the variations of E and S are intended with intensive quantities fixed. This is certainly the 
case of most thermodynamic systems also at macroscopic level (our systems included), when the volume V is 
the whole volume of the system at thermodynamic equilibrium, provided gravitational interaction or, in general, 
long range interactions can be neglected (as often assumed, at least implicitly, in thermodynamics) [l[. In this 
case the intensive quantities are intended to have values To, po and /io characterizing globally the system. 

When however gravitational effects become no longer negligible one can expect that relation ([40|) . albeit 
holding true locally (equivalence principle), be modified at a global level (with a suitable definition of the 
intensive quantities for this case), revealing thus non-extensivity. Correspondingly, relation (|4ip should be 
modified to 



S(XM) = A a S(M), 



(42) 
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where a is a constant 7^ 1 and we have indicated with M the total energy of the system. Entropy is no longer 
extensive with energy, but related to it by an Euler relation of homogeneity a =/= 1 . 

Looking at the analogy between the properties of local temperature in a lattice model with long range spin- 
spin interactions and the local temperature behaviour for gravitating systems as given by Tolman relation, the 
suggestion has been made in [l[ to define a global temperature To for gravitating systems, to be identified with 
the temperature that would be measured by an observer far away (in an asymptotically flat spacetime). 
Given this definition, thermodynamic consistency demands that for gravitating systems in self-equilibrium, or, 
in general, for closed systems not influenced mechanically by any external agent, the equation 



S = a— (43) 
Jo 

be an Euler relation of homogeneity 1/a in M, provided S here is intended as total entropy deprived of its 
component dependent on p (if present). In fact the first law (l39|) implies 1/To = dS/dM. In other terms, 
putting 7 = 1/a, wc obtain 7 — a or 



7-^, («) 

so that any 7 7^ 1 signals and measures non-extensivity. 

Black holes obey this relation, as can be verified throu gh t he explicit expressions of black hole temperature 
and entropy or, more directly, through the Smarr relation [281 ] 



(see [H, IH) where we see 7 = 1/2, so that correspondingly S(M) oc M 2 . 

Analogously (and intriguingly) for any classical gravitating system in self-equilibrium, approaching its own 
Schwarzschild radius, (conventional) entropy satisfies this same relation (|45|) giving thus again 7 = 1/2 [9(. In 
[H, equation ([44]) has been applied to a gravitating system in self-equilibrium consisting of a perfect fluid at 
constant density and the variation of 7 with M/R (being R the radius of the system) has been investigated. 

Attempting to extend the above analysis to our case, we face here the problem that, contrary to the cited 
cases, our systems are constrained to be in the assigned volume by an external pressure p ex t applied at r = R, 
so that they cannot be considered mechanically insulated. Applying however the same reasoning that leads 
to equations (|43H44)1 . we see that these equations still hold true also in case an external pressure is present, 
provided S is intended as total entropy deprived now also of the component due to a p ext 7^ 0. That is we have 



7=—, (46) 

where with So we intend the component of S obtained assuming p ext — (in addition to p = 0). 

Let us now evaluate total entropy for our systems. The Gibbs-Duhcm relation (|40|) . as applied locally, 
expressed in specific form 



p = Ts — p + pn (47) 



nvcs 



s = - (p + p-pn). (48) 

From this, considering that to deprive s of the component of p due to the presence of a p ex t ^ means to quit 
the whole p (i.e. p ex t — implies p = 0, for the fluids under study if, as at present, they are constrained to be 
in a finite volume), we have 

R R 2g-l 2g-l 
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where y = p(0) and /, h are determined by this initial condition with C = 1 (and with assigned g) and use of 
equations ([6][E]) is made, with ([6]) in the form 



b T 2g , (49) 



with b a constant. We further obtain 



S (y) = b^y^Air I r 2 ( — _!__ — ) " h^^&Ky r) 1/2 dr, (50) 



where the expression for So(y) depends only on the solutions for the initial conditions /(0) = 1; here use of the 
equation 

f(r) = J^h(V^r) (51) 

is made, relation that can readily be obtained from scaling property l|23p . Expressing the integration in (|5U|) in 
terms of the variable z = yr 2 we finally obtain 

So(z) = b^Sy-^-hn [ VI ( ]=-} M\/8^I) 1/2 dz, (52) 

Jo \fi(V8^)J 



where z = y R 2 (as in 

From equation (|48[) we see that s is given by three components, each in principle with its own radial profile. 
For our systems however these radial profiles are the same. In fact from equations ([6]) and JT|). we see that s, 
p/T and p/T terms depend only on T, with a same dependence, so that the same must hold for the pn/T term 
also. This assures the fin/T term, if present, has the same radial profile as the p/T term. Total entropy with 
all the components included, is then proportional to So so that the two entropies have a same scaling law with 
M. 

As regards global temperature To, from Tolman relation (J3J) and from l|49p and ((5J) we have 



T = T c V7(0) = b- l ^"(fi7r)- 1 ^ B C 1 ^ B , (53) 

where, as said, To is the temperature measured by a far away observer. We see that To is sensitive to the value 
of C, that on its own should depend on the mass-energy distribution of the whole system, constraining system 
included. Global temperature however is, like energy M or entropy S, a thermodynamics! parameter that, as 
such, must be characteristic of the chosen system it describes, irrespective of its environment. Moreover, if 
the constraining system is, as we suppose, in thermal equilibrium with the fluid contained in it, a same far 
away temperature Tq must coherently be measured, independently of the characteristics (density, size, ..) of 
the constraining system. From this implies that, once the asymptotic flatness of the metric is required, 
the actual value of C is fixed, irrespectively of the characteristics of the constraining system. We are thus led 
to compute To as the far away temperature evaluated as if the constraining system would be absent. From the 
continuity of the metric tensor component on the boundary of the cavity we have then 

f(R) = 1 - 2M/R, (54) 

so that from (151) 



C 



(1-2M/R)b f{R)9 : 



(55) 



where, as above, / is the solution with C = 1 and y — p(0). From this equation and from (|38|) and (|5ip we 
obtain the following expression for C 

C = 8iry ==- , (56) 
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FIG. 6: 7 versus z = p(0)R 2 for k = 1, i, 5 (highest: « = |, lowest: k = 1). 



so that Tq is 



Tq = b 2 B y2g 



[/nf^Afvte)] 1 / 2 



(57) 



From this expression for T and from the expression (|52[) for entropy S*, 7 as given by (|46[) can finally be 
written. Noting that M, given by equation ([38)) . can also be written as 



M = 2-Ky-% [ VshiVs^sy 9 dz 
Jo 

(from (I37p . once the integration variable is changed from r to z = y r 2 ), for 7 the expression 



7(2) 



M 



- WVMA^] j^p^ d2 



(58) 



is found. 



IV. SCALING LAWS AND RADIAL PROFILES 



FigElplots j(z) as given by equation (|58p. At each z value, we have S(M) oc M 1 / 7 with 7 given by the plot. 
We see that in the limit z — > (low densities and/or small radii), 7 — > 1, so that S(M) oc M (and oc R 3 , 

from ([57)) ) and extensivity obtains. This can analytically be seen from (JSSJ), as in this limit from (|2T>|) and (fTB|) 
we have /, /i ~ 1 and 7 reduces to 



7(2) ^ J ° - = 1. (59) 
J y z dz 

In the limit 2 — * 00 (very high densities and/or large radii), from ([52")) and ([55)) we get 



_j_ Jo ^ ^ 1+219 \ 9 

fo*^*- 1 dz = J + I = 9 + 1' 



so that asymptotically 7(2) = | , |, | respectively for g = 1, 2, 5 (k = 1, |, |). In this limit thus S(M) oc M 2 
for 5 = 1 and S(M) oc iW 3 with 1< /? < 2 for g > 1. Due to (J3TJ> or (J3HJ) this means also 5(i?) cc i? 2 for g = 1 
and S(R) oc i?' 3 for 5 > 1. Note that this limit can be achieved at any fixed R, if central energy density y is 
allowed to increase to sufficiently large values, and we remain far from the Schwarzschild mass limit. 
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FIG. 7: S/M 2 (arbitrary units) versus z for k = 1. 



The graph of 7 intersects the asymptote at each extremal point of ^jf of equation (|38|) (or of Fig|5|) . From 
our discussion about the stability of equilibrium configurations, this implies that the z region of stability has 
both d'y/dz < and 7(2) above its asymptote. The intersections with the asymptote also mean that for each 
configuration corresponding to a local maximum of 2M, the asymptotic behaviour of S(M) is locally recovered; 
for the g — 1 case this in particular implies that S(M) oc M 2 for each such finite z. This means that equilibrium 
configurations exist (stable, in the case of the first maximum), with y and R both finite, with exact S{M) oc M 2 
(and, due to extremality, S(R) oc R 2 ) local scaling. 

FigJT] illustrates this. We see that for z = z extr = 0.655 (the first extremal point, see FigO, S/M 2 has a local 
extremum, implying S/M 2 ~ const there. At the same time this (in addition to (|59p ) verifies the definition of 
7 as given by PS|) . 

These results confirm and, in some respects, extend what found for these systems in Within the stability 
interval in z, 7 decreases at increasing M/R, as said, and this is somehow analogous to the result obtained in 
[l[ , in the study of the relation between 7 and M/ R for a self-gravitating system consisting of a perfect fluid at 
constant density. The behaviour of conventional entropy S(R) oc R 2 confirms the statement that no horizon is 
needed to have an entropy that scales like the area Q,[9 ! ]; as our systems are far from the Schwarzschild limit 
(even when z — > 00) this statement is confirmed in a strong sense: to have an area scaling entropy, not even the 
approaching to horizon formation is required [|| . 

In terms of global temperature To, one can show that, in the limit z — - > 00, our systems have To c< 1/R [6], 
in addition to the area scaling law for entropy, so that the analogy with black hole thermodynamics is quite 
complete. This happens also for self-gravitating systems near to reach their own Schwarzschild radius Q. To 
see the To dependence on R, with reference to equation (|57[) we write TqR as 



T R=b-*Sy*S-l - - f —^' 2 / 7=7777; =v(y)C(z) (60) 



and in FigJ5]we plot £(z) for g = 1, 2, 5 (highest: g = 5, lowest: g = 1), together with the constant value 
asymptotically approached by £(z) in the g — 1 case revealing the 1/R law for Tq. From (jHO)) , together 
with equations ([32]) and (|33[) . the asymptotic behaviour of Toi? can analytically be inferred. We have TqR oc 
z°, z 1 / 4 , z 2 / 5 (corresponding to TqR oc R°, i? 1 / 2 , i? 4 / 5 ) respectively for g = 1, 2, 5 and the asymptote in FigJS] 



2V8tt 

For the g = 1 case note that every time z is extremal (for 2M/R), that happens every time the £(z) curve 
intersects the asymptote in FigJH the slope is not and thus Toi? does not scale as 1/R there, not even locally, 
notwithstanding the area scaling law for entropy. This is expected, as Tq(z) = dS(j 1 / dM and Sq(M) oc M 1 ' 7 ^^", 
so that, near each extremal z, the 1/R law can set in only if 7(2) ~ const there, that is when z — > 00. We have 
thus that only in the limit z — > 00 the thermodynamic analogy among the systems of present study and the 
black holes and the classical systems of [9| extends to both S and T scaling laws. 

Let us consider now the radial profiles of entropy. We are mainly interested in their relation with the strength 
of the gravitational interaction, as given by the parameter M/R, and in the shape they take when area scaling 
law for entropy obtains, in particular in the limiting case z —> 00 in which the 1/R scaling law for temperature 
also applies. From 
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FIG. 8: £ (representing TqR, see text) versus z for k = 1, |, | (highest: « = |, lowest: k = 1). 
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FIG. 9: Entropy local density s, normalized to its value at the centre s(0), versus radius (normalized to 1) for different 
z values (different M/R values), z = 0, z "^ tr , z ex tr- The straight line is for z = (M/R = 0); the steepest curve is for 
z — z ex tr (M/R is at its maximum). 



5' 



4irr 2 s(r)Vh dr = 



4:TTr 2 a(r) dr, 



(61) 



in addition to the entropy local density s, an entropy "radial" density a is defined (in order to compare with 
As already discussed, the three components of s as given in ([45]) have the same radial profiles. We can 
then restrict the study of entropy radial profiles to the component s = p/T alone. From (|49p and (|5ip we have 



s(r) = b^ 29 p(r)- 



i /r , 2g-l 

b 1/2g y—fi(^ 



X{y)i>z(x), (62) 



<j(r) 



x(2/)w z (x), 



(63) 



where x = ^ covers the interval [0,1] and to z {x) — i\> ' z (x)hi(\/ '8tt 'z x) 2. Being /1 monotonically increasing, we 
expect Tpz(x) to be monotonically decreasing, becoming steeper with growing z. 

For the g = 1 case, the most interesting one as regards entropy scaling laws, Figures [9l and [TOl report respec- 
tively ip z (x) and uj z {x), for z = 0, Zs g tr , z e2: tr = 0.655, where z ex t r is the first extremal point, corresponding 
to the higher edge of the stability interval in z. The straight line corresponds to z = 0, while the lowest curve 
corresponds to z = z extr - From FigO (or relation (|58|) ) we recall that M/R increases with z, from at 
z = 0, up to its maximum allowed value for z = z extr . s(r) or a(r) profiles become thus flat as the strength 
of gravitational interaction is negligible, as expected, while they become decreasing with r when M/R 7^ 0, 
being increasingly steeper with M/R growing. We have thus that the weight of the internal regions in the 
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FIG. 10: Entropy "radial" density a (see text), normalized to its value at the centre 0"(O)), versus radius (normalized to 
1) for the same z values as the previous Figure. Straight line: z — 0; lowest curve: z — z ex tr- 
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FIG. 11: r 2 s (oc entropy of the slice at radius r) versus r with s normalized to its value at the centre s(0) and r normalized 
to 1, for the same z values as previous Figures. Uppermost curve: z — 0; lowest curve: z — z extr - 



determination of total entropy grows with M/R. This can be seen also in FigfTTl where the quantity x 2 ip z (x) 
is reported (oc 4irr 2 s(r)dr, the entropy on the slice at radius r). The uppermost curve corresponds to z = 0, 
while the lowest one to z — z extr - The weight of external slices is progressively decreasing when M/R grows 
from M/R = (case with the expected quadratic dependence on r of the entropy of slices of that radius) to 
its maximum value. The behaviour of our system is at variance in this respect with the self-gravitating system 
consisting of a perfect fluid with constant density, studied in [lj. In that case the weight of external layers is 
shown to increase with M/R, as reported by Fig. 2 of [l[, the equivalent of our FigfTUl 

For z > z extr the curves for s(r) and <r(r) become increasingly steeper and it is interesting to explore what 
is the shape the radial profiles acquire in the limit z — > oo; as we know in fact in this limit for the case g = 1 we 
have both the area scaling law for entropy and the 1/R law for temperature. For z — ► oo, from (JHU) and (|35|) 
we find 



2m . 

Ui'KZjsXn 

r 



-2+i 



(64) 



and from (|32[) a same behaviour for w z {x). This means that in this limit s(r) and <r(r) behave like 1/r if g = 1 
and like l/r s with 1 < S < 2 for any other allowed g; in the limit g — > oo (dust), 5 = 2. We see thus that for our 
systems, when we are at conditions in which both the area scaling for entropy and the 1/R law for temperature 
set in, the entropy radial profile goes as 1/r, with a larger weight from internal layers in the determination of total 
entropy, as compared to the no-gravity case. This is at variance with what happens for self-gravitating classical 
systems that reach the above scaling laws for entropy and temperature while approaching their Schwarzschild 
radius, for which it is demonstrated that in this limit all entropy lies on the surface [9(. 
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V. CONCLUSIONS 



In conclusion we have seen that the area scaling law for entropy as well the 1/R law for temperature 
(the celebrated properties of black hole thermodynamics) can set in also for classical finite systems even at 
circumstances very far from the formation of an horizon. We have shown moreover that when this happens, 
entropy radial profile may be completely different from what found for the systems that gain the above scaling 
laws while approaching horizon formation, in which case entropy lies on their surface @. 

We have studied the onset of these properties for systems consisting of perfect fluids with linear equation of 
state p — up, constrained in a finite volume by an external pressure. The above properties set in at limiting 
conditions of causality [6] (k = 1) and size/energy density (R and/or p — > oo), being however safely far from 
horizon formation limit. These results have been obtained generalizing the methods introduced in [l[ for the 
study of non-extensivity of thermodynamical parameters in gravitating systems, to the case in which an external 
pressure is applied. The entropy radial profile we find in the limiting circumstances in which the mentioned 
scaling laws for entropy and temperature set in is s(r) oc 1/r, so that a heavy contribution is given by the 
internal layers, contrary to what happens for the systems [9|. This is thus a case in which the whole volume 
contributes to the total number of degrees of freedom of the system (with a heavier specific contribution from 
internal regions), being however this number scaling like the area. 
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